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Abstract. This is an investigation on the dynamical screening of an atom 
confined within a fullerene of finite width. The two surfaces of the fullerene 
lead to the presence of two surface plasmon eigenmodes. It is shown that, in the 
vicinity of these two frequencies, there is a large enhancement of the confined 
atom's photoabsorption rate. 
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1. Introduction 

Currently there is much interest in plasmon excitations and their effects on the 
properties of the object (see [1-6] and references therein). Plasmons are the collective 
excitation of electrons. They are found in strongly correlated systems such as 
nanotubes, fullerenes (not limited to the carbon variety) and metallic clusters. The 
endohedral atom, an atom confined within a fullerene, is another area of interest 
[5,7-9] and particularly metallofullerenes [10-14]. Other forms of confinement include 
confinement as dopant in metals [15] and confinement within hollow structures such 
as nanotubes and inorganic fullerenes. 

Connerade and Solov'yov studied the effect of electromagnetic screening on an 
atom confined within a fullerene [5]. An atom inside the spherical Cqq molecule was 
considered. They found that the screening is strongly dependent on the frequency, 
hence the term dynamical screening. Near the plasmon frequency of the fullerene there 
is a large amplification of the field within the fullerene, leading to an enhancement of 
the photoabsorption cross section of the endohedral atom. In their work they treated 
the fullerene as an infinitely thin spherical shell. Connerade and Solov'yov used the 
hydrodynamic formalism outlined in [6]. 

Recent work by Scully et al [16] and by Reinkoster et al [17] demonstrate the 
existence of a second giant resonance in the single photoionisation cross section of 
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Figure 1. The atom, indicated by the solid circle, is confined at the centre on 
the fullerene molecule of average radius R and thickness AR. The inner and the 
outer radii of the fullerene are given by Ri and R2 respectively. This system is 
exposed to an external field -E C xt- 



multiply ionised Cqo molecules. The first resonance is known to correspond to the 
surface plasmon of the fullerene. The nature of the second resonance has been 
discussed by Scully et al in the original paper and also by Korol and Solov'yov in 
the comment [18]. 

This is an investigation into the effects of a fullerene that has finite width, and 
in particular, focusses on the dynamical screening factor studied in [5]. 

2. Theoretical Framework 

Consider a system where an atom is confined at the centre of a spherical fullerene. This 
fullerene has an average radius of R and a thickness of AT? as depicted in Figure Q] 
When this system is exposed to an electromagnetic field, the presence of the fullerene 
dynamically screens the confined atom - it experiences a field that is either enhanced 
or suppressed depending on the field frequency to. The result is that, for the same 
external field, the photoabsorption rate of the confined atom differs from that of the 
free atom. A dynamical screening factor T = T(uS) may be defined to relate the 
photoabsorption cross sectons of these two atoms: 



In this work the problem has been studied using a dielectric shell model, first 
applied to fullerenes by Lambin et al [19] to study the two surface plasmon modes. 
The fullerene is modelled as a spherical shell of dielectric function e. It has an inner 
and an outer radius of R\ and R2 respectively, with R being their average. The atom, 
located at r = is characterised by the polarizability a a (w) and the average radius 
J? a , which is assumed to be small in comparison to R\^- This system is exposed to a 
monochromatic electromagnetic wave. Using the dipole approximation and neglecting 
the magnetic part, this is treated as an external uniform electric field E cxt (t) = E e lult . 

It is assumed that the quantities discussed here depend on the frequency w of 
the external field, i.e. all the relations are written for the Fourier components of the 
corresponding quantities. 

To calculate the screening factor one must find the absorption rate of the 

confined atom and compare it with that of the free atom. In the dipole limit, an object 
exposed to an electromagnetic wave will absorb its energy with a rate Q given by [20] 




(1) 




(2) 
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Here the integration is carried out over the volume of the considered system. E(r) is 
the electric field intensity at the position r and dD(r) is the dipole moment of the 
integration element at that position. Thus to calculate Q a , the atom's photoabsorption 
rate, it is necessary to obtain a precise expression for the field experienced by the 
confined atom. 



2.1. Calculating Electric Field at the Atom 

When the system is exposed to the external electric field, both the fullerene and the 
atom become polarized. The field, which acts on the atom may be decomposed into 
the contributions described by the following. The external field polarizes the dielectric 
shell, creating an additional field. Both fields polarize the confined atom, giving rise 
to a third field. This third field causes further polarization of the shell, which further 
polarizes the atom etc. To obtain a precise expression for the field experienced by 
the confined atom, one must consider this infinite cycle of the shell (and of the atom) 
being polarized and polarizing. 

This was achieved by splitting the problem up into two parts: the dielectric shell 
in the external field (section 12. 2|) , and a point-like dipole at the centre of the shell in 
the absence of the external field (section l2.3p . Results from these two problems were 
then combined in an iterative manner to construct the final expression for the total 
field .Etot inside the fullerene (section l2"^f)) . 



2.2. Dielectric Shell in External Field 

The problem of a solid dielectric sphere in a uniform field is well known and discussed 
in many textbooks, [20,21]. The case of a dielectric shell in a uniform field is an 
extension of that problem to a spherical layer [22,23]. Using the method detailed 
in [21], one finds the electric field to be: 



E - s ■ 



E(r) = < 
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with 



A = (2 + e)(2e + l)-2(l-e) 2 C 3 



(4) 
(5) 



and £, the ratio of the inner to outer shell radii, given by 
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2.3. Polarized Atom Inside Shell 

When a field E polarizes the atom, with polarizability a a (uj), the induced dipolc 
moment d is given by: 

d = a a (u)E. (7) 

Consider the general problem of a point-like dipole d located at the centre of the 
dielectric shell in the absence of an external field. The dipole moment polarizes the 
shell. Both fields contribute to the total field intensity over the whole space. Applying 
the same approach as in section [2721 one finds the following expression for the electric 
field: 



f _ r 2 d - 3r (dr) 

~S z 



E(r) 



',d — v- 



d — 3r (dr) 



w d 



r 2 d - 3r (dr) 



r > R 2 

Ri < r < R 2 

r < Ri, 



(8) 



where 



9e „ 6(1 -e) _ 3 (2 + e) 

A' R%A ' A ' 



3 a > 0) 



. 2 ( 2 + e)(e-l) + (l-e)(2 + e)^ 
Rf 

with A defined as in ((SJ). 



2.4- Total Field at the Endohedral Atom 

Let us now consider the problem of the atom confined at the centre of the dielectric 
shell in the presence of the external field Eq . To construct the precise expression for 
the field we combine the results from the previous sections in an iterative manner 
described below. 

Let Ei refer to the electric field at r = due to the polarization of the shell and 
di refer to the induced moment of the atom in the i th iteration. Then the iteration 
scheme is as follows. 

1 st iteration: The external field .Eo induces a dipole moment in the shell, which 
creates the additional field E\. The combined field Eq + Ei = wEq, refer to |3]), 
polarizes the atom and induces the dipole moment d\ — a&wEo. 

2 nd iteration: As discussed in section 12.31 the dipole d\ gives rise to another 
induced dipole moment of the shell. This creates the additional field E2 at r = 0, 
refer to ©: 

E 2 = wdi = [a & w)wE . (10) 
This field in turn induces the dipole moment d 2 given by: 

d 2 = a Zj E 2 = (a & w)wa & E . (11) 

i th iteration: It follows that for the iterations where i > 2, the relations for Ei 
and di are given by: 
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Ei = wdi-u di = a a Ei. (12) 

The above relations are then expressed in terms of Eq and then summed over all 
iterations to give the total field due to the polarization of the system and the total 
dipole moment e£ to t of the atom. The total electrical field at the atom -Etot, is simply 
the polarization field combined with the external field. 

E tot = FE , dtot = a a FE , (13) 

where 

F = w (l + (a a w) + (a a w) 2 + ■■■) = : - - . (14) 

1 - a a w 

The field at r = 0, given by -Etot, has been modified by the system. Equation (fT3| 
shows that the electric field, which the atom is exposed to, differs from the external 
field by a factor of F = F(u>), which depends on the frequency of the external field. 



2.5. Dynamical Screening Factor 

The absorption rate of the electromagnetic wave by the system may be divided into 
two parts: absorption due to the two surfaces of the fullerene shell and due to the 
atom. From @ and (fT5|) . the atom's absorption rate is: 

Q a = |lm(I5 t * ot (0H ot ) = |lm(a a |F| 2 |£;o| 2 ) . (15) 

The difference between this absorption rate for the confined atom and that of the free 
atom is the factor of \F\ 2 , and from (jTJ) , this factor is the dynamical screening factor 
T=\F\ 2 . 

In the approximation that the radius of the atom, R a is much smaller than that 
of the fullerene the expression for T may be simplified. The polarizability of the atom, 
a a , is roughly proportional to [21]. The term a a w is therefore negligible and (fT4]) 
becomes F ra w. Thus the screening factor is: 

(16) 



^« M 2 



9c 



Further analysis of T requires the dielectric function of the shell to be expressed 
in terms of the angular frequency of the electric field, ui. This depends on the material 
nature of the fullerene. It must be emphasised that until this point the expression for 
the screening factor is independent of the 'type of material' of the shell. Henceforth 
the fullerene shall be treated as being metallic on account of the derealization of 
the valence electrons. Following the Drude model [24] and temporarily ignoring the 
damping of the plasmon oscillations, the dielectric function e is expressed as: 

e(o; = 1 - "i. (17) 

where oj p is the volume plasmon frequency, which is described by the total number of 
delocaliscd electrons N (4 valence electrons from each carbon atom) and the volume 
of the shell V = 4ir(R% - -R?)/3: 
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It can be shown that A, see ([5]), is transformed by the Drude model into: 



(19) 



where 



6 



3-p , 



P 



p = v 7 ! + 8? 3 



(20) 



The frequencies u>\ and ll>2 are the eigenfrequencies of the surface plasmon modes of 
the fullerene [4] : the plasmons at the inner and the outer surfaces of the fullerene are 
coupled oscillators and has two normal modes. In the symmetric mode characterised 
by lo\ the charge densities of the two surfaces oscillate in phase, whereas in the 
antisymmetric mode of frequency w 2 they are in antiphase. 

The screening factor (|16p can be expressed as a sum of terms resonant in ui± and 
in UJ2- 



1 



2N 



1 



1 



iTiuj 



1T2UJ 



(21) 



In the expression above damping has been introduced via the widths of the symmetric 
and antisymmetric plasmon mode, Ti and T 2 respectively. They describe the decay 
rate from collective excitation to the incoherent sum of single-electron excitations. 
In this paper we do not attempt to calculate these widths. Various methods for 
calculating the widths are described in [25-28]. 



3. Results and Discussion 



For the purposes of calculation, the widths of the two plasmons were set as a fraction 
of the corresponding resonant frequency: Ti.2 = 7^1,2- The values of 7 used here 
are 0.1, 0.25 and 0.5. It should be noted that the screening factor is now completely 
parameterised by the the number of delocalised electrons, the outer radii and the 
quantity £, the ratio of the fullerene's inner and outer radii. 

The dynamical screening factor was studied for several spherical fullcrcncs: C20, 
Ceo, C* 2 4o, and Cgeo- The thickness of each fullerene AR was set to 1.5 A, which was 
obtained by Rudel et al [29] . The average radii of the fullerenes are given in Table [1] 
together with other relevant properties. The plasmon frequencies are calculated from 
(|20p . It must be noted that u>\ for Cqq is experimentally known to be 20 eV [16]. The 
discrepancy between the calculated values and the known experimental values may be 
explained by the fact that the model is purely classical. 

The resulting plots of the screening factor for the four fullerenes are given in 
Figure [2] The x-axis is the photon energy scaled with respect to wi to aid comparison 
between the different fullerenes. It should be noted that C20 has the largest u>%. 

All the plots of the screening factor show that in the static limit the atom is 
completely shielded from the external field. At every high frequencies the fullerene 
is transparent to the light and there is no enhancement. This is as expected for a 
metallic fullerene. The profile of the screening factor for the different fullcrcncs are 
very similar and show a strong dependence on the value of 7. There is a peak in the 
vicinity of the first resonant frequency and another feature near the second resonant 
frequency due to the symmetric and the antisymmetric plasmon modes, respectively, 
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Figure 2. Plots of the screening factor for the 4 fullerenes. For each fullerene, 
the screening factor was plotted for the same set of 7. The x-axis is the photon 
energy scaled with respect to that of the first resonance. The pairs of thin vertical 
lines indicate the position of the two resonant frequencies for each fullerene. 



Table 1. Properties of the 4 fullerenes are given in the table: the average radius, 
the ratio of the inner and outer radii, and the resonant frequencies calculated 
from J20t . 



Fullerene 


R(A) 




wi (eV) 


W2 (eV) 


C20 


2.0 [30] 


0.45 


19.8 


31.7 


Ceo 


3.5 [29] 


0.65 


16.9 


33.5 


C240 


7.1 [31] 


0.81 


12.8 


35.0 


C96O 


13.8 [31] 


0.89 


9.9 


37.1 



as demonstrated in Figure [3] Their characteristics are determined by the damping 
of the plasmons, which are determined by the value of 7. For each fullerene, there is 
a well defined peak at each resonant frequency for 7 = 0.1. For larger values of 7, 
the individual contribution of the second plasmon mode is strongly suppressed by the 
damping. Here the interference terms of the screening factor become more significant 
and the second resonance is reduced to an extension of the first peak, as shown in 
Figure [3] 

From Figure [5]it can be seen that the relative heights of the first peaks are almost 
the same for the different fullerenes. In this region the contribution from the second 
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Figure 3. A plot of the screening factor for C20 (left panel) and Ceo (right 
panel) with 7 = 0.25. The individual contributions of the symmetric and the 
antisymmetric plasmon modes are also shown. 



plasmon mode is effectively constant and is small. Using the approximation that the 
maximum occurs at the first resonant frequency, one arrives at this result: 




where p is a function of £, which is fixed by the particular fullerene. This relation 
explains the large differences in J 7 max , seen in Figure [H for the different values of 7 
and also between the different fullercncs. 



3.1. Infinitely Thin Fullerene 

In the limit of the infinitely thin fullerene of radius R, £ — > 1. It is useful to note that 
(1 — £ 3 ) = (p+ 3)(p — 3)/8. Following on from this, one shows that the two resonant 
frequencies become 



lim ui\ = lim — ( 3 

£-»i £->i 6 v 



lim oj, = lim — | 3 

€-1 €-1 6 



27V 
3~R2 



(23) 



P 



Physically, the two surfaces approach each other. At the limit, the two surfaces 
merge and only the symmetric mode is possible. Here, L02 — > 00, removing the second 
resonant term, and u)\ becomes Wf, the resonant frequency for the surface plasmon of 
the infinitely thin fullerene. Thus the screening factor becomes: 



lim T = 



1 



Tito 2 



(w 2 



,2^2 



) 2 + 



(24) 



This is the result given by Connerade and Solov'yov in [5] . A comparison of the two 
screening factors is given in Figured] For Cqo loj is 20.2 eV. 

The distinguishing feature of the screening factor of the thick fullerene is the 
presence of the second peak. Recent studies [16, 17] have demonstrated the existence 
of a second plasmon peak in the photoabsorption cross section of the fullerene. 



Dynamical Screening of an Atom Confined within a Finite- Width Fullerene 



9 




Photon energy (eV) 



Figure 4. A plot comparing the screening factors for the infinitely thin fullerene 
[5] and that of the fullerene of finite thickness. This is done for Cgo with 7 = 0.25. 

4. Conclusion 

In this paper we have generalised the dynamical screening factor obtained by 
Connerade and Solov'yov [5] for an atom confined within an infinitely thin fullerene 
to that of an atom confined by a fullerene of finite width. We have taken a classical 
dielectrical approach to study the screening of an atom at the centre of the fullerene. 
We have focussed on the case where the polarisability of the confined atom is 
sufficiently small so as to not have an appreciable effect on the field inside the fullerene. 

With this classical approach we give a qualitative description of the dynamical 
screening factor. We have shown that there are two surface plasmon eigenomdes 
of this thick fullerene. They result in two frequency ranges where there are large 
enhancements of the confined atom's photoabsorption cross section. This is in contrast 
to the single frequency range for the infinitely thin fullerene. For cases of minor 
damping of the plasmons, there are two distinct peaks of the screening factor in the 
vicinity of the plasmon resonant frequencies. When damping is large, the second peak 
is suppressed and becomes an extension of the first peak. 

Despite the simplicity of this model, it provides an instructive description of the 
phenomena discussed here. This model may be further applied to study the dynamic 
responses of more complicated situations various processes. 
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